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I. INTRODUCTION 

Aging behaviour observed in systems with slow 
degrees of freedom is one of the most intriguing aspects 
in nonequilibrium physics, see Refs. \WM for recent 
reviews. This behaviour is due to relaxation processes 
which depend on the thermal history of the studied 
sample. Aging phenomena have been studied exten- 
sively in disordered systems like spin glasses or glassy 
systems, but they are also encountered in the simpler 
ferromagnetic systems In the last years remark- 
able progress has been achieved in our understanding 
of the physics of aging phenomena taking place in 
ferromagnets S- 6 - 7 - 8 - 9 - 10 - 11 - 12 - 13 - 14 - 15 - 16 - 17 - 18 - 19 - 20 - 21 - 22 - 23 - 24 - 25 

In the typical scenario the ferromagnetic system is pre- 
pared in a fully disordered state at high temperatures 
and then quenched down to temperatures T < T C1 where 
T c > is the critical temperature. When the final tem- 
perature is smaller than T c , coarsening takes place, lead- 
ing to the formation of domains with a time-dependent 
typical length. It is the slow motion of the domain bound- 
aries which is responsible for the aging processes. The 
situation is different for a quench to the critical point, 
where one is dealing with nonequilibrium critical dynam- 
ics, as ordered domains do not exist in that case. The 
typical time-dependent length is then given by the dy- 
namical correlation length £(t) ~ i 1 / 2 , where z is called 
the dynamical critical exponent. 

In recent years investigations concentrated on two-time 
quantities like dynamical correlation and response func- 
tions, as they more fully reveal the aging processes. Ex- 
amples are the spin-spin autocorrelation function C{t, s) 
and the conjugate response function R(t, s) defined by 



dynamical scaling regime, t t mlcl:o , s ^> t m i cro and 
r = t — s 3> tmicro, where t m i cro is some microscopic 
time, these quantities are expected to have the following 
behaviour 1-2 : 



C(t,s) 



b fc(t/s) , R(t,s)=s- 1 - a f R (t/ S ). (2) 



In the limit y = t/s — > oo the scaling functions fc,n(y) 
present a simple power-law behaviour 



fc(y)~y- Xc/z , f R (y)~y- Xn/z 



(3) 



with the autocorrelation-^ and autoresponse*^ exponents 
Ac and A^. In the case of a quench to the critical point 
—of a ferromagnetic system with initial short-range corre- 
lations we have a = b = 2x/z and Ac = A_r = d — Xi + x 
where d is the number of space dimensions, whereas x 
resp. Xi denotes the scaling dimension of the order pa- 
rameter resp. of the initial magnetization^ A generaliza- 
tion of dynamical scaling to local space-time dependent 
scaling yields the predictionii*i4 



C(t,s) = ( ( f>(t)0( s )) , R(t,s) 
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Here 4>(t) is the time-dependent order parameter, 
whereas h(s) is the magnetic field conjugate to <j>. The 
time t elapsed since the quench is called the observa- 
tion time and s is referred to as the waiting time. Of 
course, R(t, s) — for t < s, due to causality. In the 



fn{y) = r y 



(y-iy 



(4) 



for the scaling function of the response function R(t,s). 
This prediction was found to be in excellent agree- 
ment with Monte Carlo simulations of two- and three- 
dimensional Ising models both at and below T c m^i^^ 
Similar agreement has recently also been observed for the 
three-dimensional XY models On the other hand re- 
cent renormalization group calculations^ yielded at crit- 
ical points a correction at two loops to the prediction Q 
coming from the theory of local scale invariance^ The 
origin of this discrepancy is still not elucidated, but it 
must be noted that the simulations are based on a mas- 
ter equation, whereas the field theoretical calculations 
start from a Langevin equation. 

Recently, the fluctuation-dissipation ratio 



X(t,s) = TR(t,s)/ 



dC(t,s) 
ds 



(5) 



has attracted much attention in nonequilibrium 
systems^ At equilibrium, both C(t,s) and R(t,s) 
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depend only on the time difference t — s, and X(t, s) = I, 
thus recovering the well-known fluctuation-dissipation 
theorem. However, this is not the case far from 
equilibrium where X(t,s) =/= 1. This property led to 
numerous studies of X(t, s) in various systems with the 
aim to possibly quantify the distance to equilibrium. 
In the context of glassy systems the limit value X°° 
has been interpreted as giving rise to an effective 
temperature T e ft = T/X 00 ^ Godreche and Luck 
conjectured that in critical systems the asymptotic value 

X°° = lim I lim X(t, s)) of this nonequilibrium 

s >oo \t >oo / 

quantity characterizes the different universality classes* 9 " 
This has indeed been confirmed in exactly solvable 
cases; 8 ! 10 ' 30 in various Monte Carlo studies 1 ! 9 ! 16 ! 17 ! 25 ! 31 
as well as in field-theoretical calculations i 13 

Surfaces are unavoidable in real materials and play an 
outstanding role in nanotechnologie. It is therefore of 
importance to understand the local behaviour close to a 
surface in nonequilibrium systems. In this work we take 
a step in that direction and extend the investigation of 
aging phenomena in critical systems to the semi-infinite 
geometry, focusing on aging processes taking place in the 
vicinity of the surface. 

It is well known that the presence of symmetry- 
breaking surfaces changes local quantities, leading to sur- 
face critical behaviour which differs from that of the bulk 
(see Refs. I32H33I for recent reviews). For example, the 
static critical surface pair correlation function behaves 
as \p — p'\ 2x1 , where for a <i-dimensional model p is a 
(d — 1 )-dimensional vector parallel to the surface. The 
scaling dimension x\ of the surface order parameter dif- 
fers in general from the scaling dimension x of the bulk 
order parameter. Of further importance is the fact that, 
depending on the surface interactions, different surface 
universality classes may be realized for a given bulk uni- 
versality class. An important case is that of the three- 
dimensional semi-infinite Ising model where three differ- 
ent surface universality classes are encountered at the 
bulk critical temperature. 

The present work studies the effect of surfaces on ag- 
ing processes with a non-conserved order parameter in 
the dynamical scaling regime t,s^$> t m i cro , t — s ^> 
tmicro- We thereby discuss not only two- and three- 
dimensional semi-infinite Ising models, but also ana- 
lyze aging phenomena taking place in the Hilhorst-van 
Leeuwen modeli^i^ 5 . The latter model is a semi-infinite 
two-dimensional Ising model with an extended surface 
defect. In this model the surface scaling dimension Xi 
has been shown to vary continuously as a function of the 
defect amplitude in the interesting case that the surface 
defect is a marginal perturbation. This makes it possi- 
ble to systematically study at criticality the impact of 
surfaces on local two-time correlation and response func- 
tions as well as on related quantities as local fluctuation- 
dissipation ratios. 

This work is also meant to close a gap in the study of 
critical dynamics at surfaces. Prior works either treated 



equilibrium critical dynamics near surfaces^^LiS or fo- 
cused mainly on the effect a surface has on nonequilib- 
rium critical dynamics immediately after the quench to 
the critical pointi^ 9 ^ ^^ 2 ^ 3 " To our knowledge the dy- 
namical scaling regime has not yet been investigated in 
semi-infinite critical systems. 

The paper is organized in the following way. In Sec- 
tion II the different semi-infinite models are introduced 
and some important facts on equilibrium and dynamic 
surface critical behaviour are reviewed. Section III deals 
with surface two-time autocorrelation functions in the 
dynamical scaling (i.e. aging) regime, whereas Section IV 
is devoted to the corresponding study of the thermore- 
manent magnetization. These quantities are examined in 
the Hilhorst-van Leeuwen model for various values of the 
defect amplitude. In the three-dimensional semi-infinite 
Ising model, both the ordinary transition (where the bulk 
alone is critical) and the special transition point (where 
bulk and surface are both critical) are considered. In 
Section IV numerically determined scaling functions de- 
rived from the intergrated responses are compared with 
predicted scaling functions coming from the theory of lo- 
cal scale invariance. Section V is devoted to the analysis 
of surface fluctuation-dissipation ratios. It is shown that 
the limit value of this ratio depends on the value of the 
surface scaling dimension x\. Interpreting this asymp- 
totic value as giving rise to an effective temperature, it 
follows that for a given bulk universality class different 
surface effective temperatures may be obtained by vary- 
ing the value of some local model parameters. Finally, 
Section VI gives our conclusions as well as an outlook on 
open problems. 



II. CRITICAL BEHAVIOUR NEAR SURFACES 

A. Models and surface phase diagrams 

All the models studied in this work are defined on 
e?-dimensional semi-infinite (hyper)cubic lattices. Every 
lattice point is characterized by an Ising spin which takes 
on the values ±1. For the pure, perfect Ising model in 
absence of external fields, the Hamiltonian is given by 

H = ai<T o ~ Js a i a r ( 6 ) 

bulk surface 

The first sum is over nearest neighbour pairs where at 
most one spin is a surface spin, whereas the second sum 
is over nearest neighbour pairs with both spins lying in 
the surface layer. Jt resp. J s is the strength of the bulk 
resp. surface couplings. All couplings are supposed to be 
ferromagnetic, i.e. Jb, J s > 0. 

The surface phase diagrams of the two- and three- 
dimensional semi-infinite Ising models © are well 
establishedj 3 ^ 3 - In two dimensions only one surface uni- 
versality class, the so-called ordinary transition, is en- 
countered at the bulk critical point for all values of the 
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surface couplings. This transition is characterized by the 
value X\ = 1/2 of the scaling dimension of the surface 
order parameter. Here and in the following I follow the 
usual convention and characterize surface quantities by 
the index 1. The surface phase diagram is more interest- 
ing for the three-dimensional model. If the ratio of the 
surface coupling J s to the bulk coupling Jt>, r — J s /J{,, 
is sufficiently small, the system undergoes at the bulk 
critical temperature T c an ordinary transition, with the 
bulk and surface orderings occuring at the same temper- 
ature. Beyond a critical ratio, r > r sp w 1.50 for the 
semi-infinite Ising model on the simple cubic lattice^^ 
the surface orders at the so-called surface transition at a 
temperature T s > T c , followed by the extraordinary tran- 
sition at T c . At the critical ratio r sp , one encounters the 
multicritical special transition point, with critical surface 
properties deviating from those at the ordinary transition 
and those at the surface transition. The present work 
exclusively treats aging processes taking place at the or- 
dinary transition and at the special transition point. For 
these cases the scaling dimension x\ takes on the values 
1.26 (ordinary transition) and 0.376 (special transition 
point). 

The Hilhorst-van Leeuwen model^S is a critical two- 
dimensional semi-infinite Ising model with an extended 
surface defect due to inhomogeneous couplings. Con- 
sidering a square lattice, one has couplings with a con- 
stant strength J\ in the direction parallel to the surface, 
whereas the strength of the couplings varies perpendicu- 
lar to the surface as a function of the distance / to the 
surface: 

J 2 (0-J 2 (oo) = ^. (7) 

with A = AT c sinh(2J 2 (cx))/r c )/4. 46 In this work J x is 
set equal to J 2 (oo). This extended perturbation is irrel- 
evant for w > 1, yielding the same critical behaviour as 
the homogeneous semi-infinite system. For oj < 1 and 
A > 0, however, the perturbation is relevant and a spon- 
taneous surface magnetisation is observed at the bulk 
critical point. In the present context the most interest- 
ing case is obtained for u> = 1, where the perturbation 
is marginal for A < 1 (for A > 1 one again observes 
a spontaneous surface magnetisation). Indeed, exact re- 
sults show that the scaling dimension x\ is then a con- 
tinuous function of the parameter A with 

xt = \{l-A) (8) 

and A < 1. This intriguing behaviour, which has at- 
tracted much interest in the past, gives us the possi- 
bility to continuously change the surface critical expo- 
nents for a given bulk universality class (that of the two- 
dimensional Ising model), thereby making a systematic 
study of aging processes at critical surfaces possible. 



TABLE I: Literature values of bulk and surface scaling dimen- 
sions x and xi as well as of the nonequilibrium dynamical 
bulk exponents Xi and z, for both the two- and the three- 
dimensional Ising models. OT: ordinary transition, SP: spe- 
cial transition point. 





X 


21 


Xi 


z 


d = 2 


1/8 


1/2 


0.53 


2.17 


d = 3, OT 


0.516 


1.26 


0.74 


2.04 


d = 3, SP 


0.516 


0.376 


0.74 


2.04 



B. Critical dynamics near surfaces 

Besides changing the local static critical behaviour, 
surfaces have also an effect on critical dynamics. This has 
been studied in semi- infinite extensions of the well-known 
bulk stochastic models^ as for example model A (purely 
relaxational without any conserved quantities^S^i or 
model B (with conserved order parameter) A central 
aspect of the works on dynamic surface critical behaviour 
concerns the possible classification of the distinct surface 
dynamic universality classes. 

It is important to note that no genuine dynamic sur- 
face exponent exists i^SiSi All exponents describing the 
surface critical behaviour of dynamic quantities can be 
expressed entirely in terms of static bulk and surface ex- 
ponents and the dynamic bulk exponents z and X{. The 
values of the critical exponents that are of interest in 
the following are listed in Table [I] for the two- and the 
three-dimensional Ising models. Consider as an exam- 
ple the dynamic surface spin-spin autocorrelation func- 
tion. At equilibrium it decays for long times as t~ 2x± / z , 36 
whereas out-of-equilibrium one expects from scaling ar- 
guments and field-theoretical calculations the power-law 
behaviour t~ Xl ^ z with the surface autocorrelation expo- 
nent 

Xi=X c + 2( Xl - x) (9) 

and t > 1 being the time elapsed since the quenchiSLiSi 
This prediction has up to now only been verified in semi- 
infinite Ising models. 

Out-of-equilibrium studies of dynamic surface crit- 
ical behaviour are very scarcei^^S^^ 2 ^ In Ref. |3^ 
short-time critical relaxation was analysed at surfaces 
by preparing the system at high temperatures with a 
small initial surface magnetization. As in the bulk case, 
a power-law behaviour of the surface magnetization is 
observed at early times, governed by the nonequilibrium 
exponent 9\ = (xi — x\)jz. Recently the two-time sur- 
face autocorrelation function Ci(t, s) was analysed in the 
short-time regime t — s <C s. 43 Using scaling arguments 
it was predicted that in the case Xi < x\ a new effect, 
called cluster dissolution, takes place, which should lead 
to an unconventional, stretched exponential dependence 
of the short-time autocorrelation: 

Ci (t - s) ~ exp(-C(t -s) K ) (10) 
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with 



It follows that for t — s <C s a stationary behaviour 
should be observed, as the autocorrelation in this case 
only depends on the time difference t — s. These predic- 
tions were verified by simulating the dynamical evolution 
of the Hilhorst-van Leeuwen model for various values of 
x\. Interestingly, this stretched exponential behaviour is 
also observed in the three-dimensional semi-infinite Ising 
model at the ordinary transition— where Xi = 0.74 and 
xi = 1.26. For the case Xi > Xi, on the other hand, 
no stationary behaviour is predicted for t — s <C s. It is 
worth noting that dynamical surface response functions 
have not yet been analyzed in critical semi-infinite sys- 
tems. 

In the dynamical scaling regime (t,s,t — s t m icro) 
it is expected that the surface autocorrelation function 
Ci(t, s) and the surface autoresponse function R\(t, s) 
exhibit a behaviour similar to that of the corresponding 
bulk quantities, but with surface exponents replacing the 
bulk exponents: 

d(t,s) = S - hl / Cl (i/ S ) (12) 
Ri(t,s) = s- 1 -* 1 f Rl (t/a) (13) 

The scaling functions fd (y) and f Rl (y) display again a 
power-law behaviour in the limit y — ► oo: 

fcM-v-* '' > fRi(v)~v~ xy '- (i4) 

In the case of a quench to a critical point from an initially 
uncorrelated state, one expects 

= = Ai (15) 

where Ai is given by Eq. @. Furthermore, one should 
have tha1i2& 

ai = bi = 2xi/z (16) 
where x\ is again the surface scaling dimension. 



III. SURFACE AUTOCORRELATION IN THE 
AGING REGIME 

The scaling behaviour (|12|) of the surface autocorrela- 
tion function is predicted to be valid in the dynamical 
scaling regime for all values of x\. In the short-time 
regime, however, one has to distinguish whether x% > x%, 
so that the usual domain growth mechanism prevails (this 
is the same situation as encountered in bulk systems), or 
whether Xi < Xi, in which case cluster dissolution takes 
placed leading to the stationary stretched exponential 
behaviour (|1(J|) . In the domain growth regime, the scaling 
form HI 2(1 is expected to be valid even for t — s < s. 



These predictions are confronted in the following with 
results obtained in extensive Monte Carlo simulations 
of the Hilhorst-van Leeuwen model and of the three- 
dimensional semi-infinite Ising model. For the former 
model, various values of the amplitude A were consid- 
ered, with the value of the scaling dimension of the sur- 
face magnetization ranging from Xi = 1/4 to x % = 1, 
see Eqs. J7J and ©■ The homogeneous two-dimensional 
semi-infinite Ising model is recovered in the special case 
A = and x\ = 1/2. The different systems were pre- 
pared at t = at infinite temperature (i.e. completely 
uncorrelated initial state) and then quenched down to the 
critical point. Applying heat-bath dynamics, the dynam- 
ical evolution of the system was studied for times t > 0. 
Heat-bath dynamics was chosen so that no macroscopic 
quantities were conserved. After the waiting time s the 
time dependence of the dynamical surface autocorrela- 
tion function Ci(i, s) was measured, with 

<m*.*) = ^ E fa(*w«)> ( 17 ) 

i£surface 

and t > s. The sum in ((17(1 is over all surface spins, 
N is the total number of surface sites, whereas <7j(t) 
is the value of the spin located at the surface site i 
at time t. The brackets indicate an average over ther- 
mal noise. Typically, the average has been taken over 
more than 10000 different realizations with different ran- 
dom numbers. In the actual simulations, systems with 
300 x 300 spins were simulated in two dimensions, with 
periodic boundary condition in the direction parallel to 
the surfaces and free boundary condition in the non- 
homogeneous direction. For the three-dimensional Ising 
model, systems with typically 40 layers with 60 x 60 spins 
per layer were studied. Here, periodic boundary condi- 
tions were used in both directions parallel to the sur- 
faces. In both cases some simulations for other system 
sizes were also done, in order to check against finite-size 
effects. All the data discussed in this work are free from 
this kind of effects. 

Figure^gives a first impression of the impact that the 
surface scaling dimension x\ has on the surface two-time 
autocorrelation function. The Figure displays autocor- 
relation functions obtained in the Hilhorst-van Leeuwen 
model for various values of x\, ranging from x\ = 1/4 
(top) to x\ — 1 (bottom), with the waiting time fixed at 
s = 25. The grey dashed line is the function obtained 
for A = 0, i.e. for the homogeneous two-dimensional 
semi-infinite Ising model. The data show that when the 
value of x\ increases, the decay of the correlation be- 
comes faster. This is due to two different effects: (a) in 
the long-time limit, the scaling function fci(t/s), which 
is proportional to C\(t, s), displays a power-law decay 
with the exponent Xi/z, where Ai increases with x\ (see 
Eq. J??}), and (b) at short times with X\ > Xi, the initial 
behaviour is that of a stretched exponential where the 
exponent k again increases with x\ j*» see Eq. (II X I) . 

The scaling function /ci(V s ) = s 2xi / z C\(t, s) is 
shown in Figure [3 for three different values of the seal- 
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FIG. 1: Dynamical surface spin-spin autocorrelation functions 
vs t — s for s fixed at 25, as obtained in the Hilhorst-van 
Leeuwen model. The different curves correspond to different 
values of the surface scaling dimension: xi = 1/4, 3/8, 1/2, 
5/8, 3/4, 7/8, 1 (from top to bottom). The dashed grey line 
is obtained for the pure homogeneous two-dimensional Ising 
model with X\ = 1/2. 



TABLE II: Comparison of numerically determined values of 
Ai/z with the theoretical expectations \\j z has been de- 
termined from both the autocorrelation (C) and the autore- 
sponse (R) functions. The errors result from averaging over 
different waiting times. The available numerical data do not 
permit to reliably determine Xi/z for x\ > 3/4 in the case of 
the Hilhorst-van Leeuwen model. HvL: Hilhorst-van Leeuwen 
model, OT: ordinary transition, SP: special transition point. 



HvL 


A 


Xl 


expected 


\i/z(C) 


Xi/z(R) 




0.50 


1/4 


0.85 


0.84(1) 


0.84(1) 




0.25 


3/8 


0.96 


0.96(1) 


0.95(1) 







1/2 


1.08 


1.09(1) 


1.12(2) 




-0.25 


5/8 


1.19 


1.22(2) 


1.21(2) 




-0.50 


3/4 


1.31 


1.30(2) 


1.35(3) 


d = 3 


OT 


1.26 


2.10 


2.10(1) 


2.18(3) 


d = 3 


SP 


0.376 


1.22 


1.16(2) 


1.24(2) 



ing dimension x\. In all cases a perfect data collapse 
is achieved when plotting fci{t/s) for different waiting 
times. As x\ > Xi for x\ =3/4, one has in that case a sta- 
tionary stretched exponential behaviour at early times^ 
and therefore the scaling is not expected to work at short 
times. This is indeed the case, but this effect is not visi- 
ble on the scale used in Figure El (see, however, the inset 
in Figure 21 below). The range of waiting times shown 
here is comparable to that of Ref. 0, where aging pro- 
cesses taking place in critical bulk systems were investi- 
gated. The dashed lines illustrate the expected long-time 
behaviour l|14|) of the scaling function, which is in com- 
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FIG. 2: Logarithm of the scaling function fci(t/s) = 
g 2xi/z (^(t, s) as function of the logarithm of the scaling vari- 
able t/s. The data shown have been obtained in the Hilhorst- 
van Leeuwen model for three different values of the surface 
scaling dimension x\ and various waiting times s. The dashed 
lines indicate the expected power-law behaviour for t/s S> 1. 
Here and in the following, errors are comparable to the scat- 
tering of the data. 



plete agreement with the numerical data, see Table [H] 
Unfortunately, no theoretical prediction of the complete 
functional form of the scaling function exists. There has 
been some recent progress in the derivation of an analyt- 
ical expression for fc in the bulk system^S^^ but these 
achievements are up to now only applicable to the spe- 
cial case z = 2 not realized in the critical systems we are 
interested in. 

In the three-dimensional semi-infinite Ising model dif- 
ferent surface universality classes are encountered when 
changing the strength J s of the surface couplings. This 
is not the case in two dimensions, where one has only 
the ordinary transition for all values of J s . Neverthe- 
less, the dynamical two-time autocorrelation function is 
affected by a change of the value of J s . This is illustrated 
in Figure |3| where the logarithm of the scaling function 
is plotted for three different values of r = J s /Jb- 5/4, 
1, 3/4, and two different waiting times. One remarks 
that C\ (t, s) decreases with increasing value of r for a 
given value of t/s. However, when multiplying Ci(t, s) 



6 




FIG. 3: Logarithm of the scaling function fd(t/s) — 
s 2xi ' z Ci(t, s) obtained in the two-dimensional semi-infinite 
Ising model as function of ln(f/s) for three different values of 
the strength J s — r Jt of the surface couplings, with r = 5/4 
(top), 1, and 3/4. The black (grey) lines correspond to wait- 
ing times s = 100 (50). 
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by a constant so that the different curves coincide for 
large values of t/s, one observes that the rescaled func- 
tions are identical in the whole aging regime. It is only 
at short times with t — s < s that slight deviations are 
observed. This is not surprising, as the short-time be- 
haviour of the dynamical autocorrelation function is not 
expected to be universal in the present case. Concerning 
the aging regime, one may conclude that the scaling law 
(|12H should in fact read 

c 1 (t, s ) = n{j s ) s - 2 ^' z f Cl (t/ s ) (is) 

where only the amplitude Q depends on the strength J s 
of the surface couplings. The same conclusion holds for 
values of A ^ 0. 

Finally, Figure0]shows the scaling functions fc x deter- 
mined in the semi-infinite three-dimensional Ising model 
for two different surface universality classes: the ordinary 
transition (with J s = Jf,) and the special transition point 
(with J s = 1.5 Jb). In both cases data collapse is achieved 
when inserting the correct values of the surface scaling 
dimension, see Table [I] The inset illustrates the short- 
time behaviour. Recall that at the ordinary transition 
x\ = 1.26 > Xi, whereas at the special transition point 
xi = 0.376 < Xi. Following the scaling arguments given 
above, a stationary behaviour is expected in the former 
case, yielding at short times deviations from the scaling 
behaviour (|12f> . whereas in the latter case the same scal- 
ing behaviour should also prevail at short times. This is 
indeed what is observed. 



FIG. 4: The same as Figure |3 but now for the three- 
dimensional Ising model at the ordinary transition (OT) and 
at the special transition point (SP). The dashed lines indi- 
cate again the expected power- law behaviour for t/s 3> 1. 
Inset: short-time behaviour showing the predicted deviations 
from scaling at the ordinary transition, whereas at the special 
transition point the scaling behaviour 1121 is also observed at 
short times. 



IV. SURFACE THERMOREMANENT 
MAGNETIZATION 

In order to compute the surface thermoremanent mag- 
netization we follow Barrat 6 and prepare the system in an 
uncorrelated initial state before quenching it down to the 
critical point in presence of a small binary random field 
hi = ±/i. After the waiting time s the field is switched 
off and the surface (staggered) magnetization 

Mi(t,s) = ± ( h *Mt))/Tc (19) 

iGsurface 

is measured, with t > s being the time elapsed since 
the quench. Here (■ ■ •) indicates again an average over 
the thermal noise whereas the bar means an average over 
the random field distribution. The sum in Eq. (|19fl is 
restricted to the N surface sites. 

It is very expensive to obtain good data for the surface 
integrated response function at larger values of the ratio 
t/s. This is due to the value of the surface scaling dimen- 
sion, which governs to a large extent the behaviour in the 
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aging regime and which in general considerably exceeds 
the value of the scaling dimension in the bulk. In the 
simulations I went up to observation times t = 21 s and 
typically averaged over 500000 different realizations for 
every waiting time considered. The sizes of the systems 
are the same as those used for the study of the surface 
autocorrelation function. 

The local thermoremanent magnetization is related to 
the surface response function by 



Mi0,s) =h duRi(t,u) 



(20) 



Mi (t, s) should have the following scaling behaviour in 
the dynamical scaling regime 



M 1 (t,s)/h = s- 2x ^f Ml (t/s) 



(21) 



with the scaling function /mi> Expression (|21|l results 
from inserting the expected behaviour i|13|) of the surface 
response function into Eq. (|20|l . 

Exact expressions for the scaling functions and /mi 
are obtained by assuming that local scale invarianceii* 1 ^ 
holds. As shown in the Appendix, one then obtains for 
the surface autoresponse the expression 



and, after integration, 

htM = r y- x ^ z 
2 -Fi(l + 2xi/z, \i/z - 2x\/z\ Xi/z 



1) 



-1— 2xi/z 



(22) 



2xi/z + l; 2 T 1 )(23) 



where 2F1 is a hypergeometric function. Expressions Ill'2ll 
and H23JI are similar to those obtained for bulk systems, 14 
but now the surface scaling dimension x\ and the surface 
autoresponse exponent Ai replace the corresponding bulk 
quantities. 

The scaling form 1)21(1 results by naively inserting Eq. 
(|13|l into Eq. 120|) without taking into account that the 
condition for the validity of (|13|l is violated at the lower 
integration bound. It was shown iniS that in case the 
system undergoes phase-ordering at temperatures below 
T c the leading correction-to-scaling term can usually not 
be neglected. However, at a critical point this correction 
term, which is of the form s~ Xl ^ z gu x (t/s)r& where <?Mi 
is a scaling function known in the framework of the theory 
of local scale invariance, should not have a sizeable effect, 
as X±/z is much larger than 2x\/z. (In case of a bulk 
system, the local exponents Ai and x\ have of course to 
be replaced by the corresponding bulk exponents.) This 
has been confirmed in the study of aging phenomena in 
critical bulk systems^ 1 - and it is also the case for the 
surface thermoremanent magnetization discussed in the 
following. 

The behaviour of the autoresponse function in the 
short-time regime t — s < s has not yet been discussed 
in the literature. Some conclusions regarding the autore- 
sponse function may, however, be drawn by remarking 
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FIG. 5: The thermoremanent magnetization as function of 
the time difference t — s in the Hilhorst-van Leeuwen model 
with xi = 1 and h = 0.1. A stationary regime is observed at 
early times as M\{t, s) only depends on t — s. In this figure 
error bars are comparable to the symbol sizes. 



that for t — s <C s the system is in local equilibrium, 
yielding a fluctuation-dissipation ratio X(t, s), see (0, 
that takes on its equilibrium value l. 3 This ratio shows 
strong deviations from its equilibrium behaviour only at 
later times when leaving the quasiequilibrium regime, as 
discussed in the next Section. It may therefore be con- 
cluded that the autoresponse function should exhibit at 
short times a stationary behaviour every time the au- 
tocorrelation function displays this kind of behaviour. 
This is indeed the case, as illustrated in Figure [5] for the 
Hilhorst-van Leeuwen model with X\ = 1. 

The main findings of this Section are summarized in 
Figures and Q showing the rescaled thermoremanent 
magnetization in the Hilhorst-van Leeuwen model for dif- 
ferent values of the surface scaling dimension as well as in 
the three-dimensional Ising model at the ordinary transi- 
tion and at the special transition point. One remarks that 
in all cases the thermoremanent magnetization shows a 
perfect scaling behaviour and that the scaling function 
/mi(£, s) — s 2xi / z Mi(t, s) for t — s > s indeed only de- 
pends on the scaling variable t/s. Deviations from scaling 
are observed for t — s < s when Xi < x%, as expected. 
The scattering of the data in these figures being a mea- 
sure of the error bars, it is obvious that reliable data are 
more difficult to obtain for larger values of x\. This fact 
precludes the study of the autoresponse function in the 
aging regime for values of x\ much larger than 1. 

The full lines in Figures El and are obtained from 
the theoretical prediction l(23|l by inserting the numerical 
values of the scaling dimensions x\ and Xi. The non- 
universal amplitude tq has thereby be fixed by requiring 
that the analytical curve and the numerical data coincide 
for large values of t/s. In cases where Xi > X\ the theo- 
retical curves perfectly describe the numerical data over 
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FIG. 6: Logarithm of the scaling function of the thermorema- 
nent magnetization as function of ln(t/s), as obtained in the 
Hilhorst-van Leeuwen model for three different values of the 
surface scaling dimension xi and for different waiting times. 
The strength of the binary random field was h — 0.1. The 
full grey lines follow from the theoretical prediction 12311 . 
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FIG. 7: The scaling function of the thermoremanent mag- 
netization vs t/s, as obtained in the three-dimensional semi- 
infinite Ising model at the ordinary transition (OT) with J s = 
Jb and at the special transition point (SP), with h = 0.05. 
The full grey lines are the theoretical predictions lY2>>t coming 
from the theory of local scale invariance. 



the whole range of values of the ratio t/s. Deviations 
are observed in the short-time regime t — s < s when 
Xi < x\. This is obvious in Figure for x\ — 3/4, but it 
is also the case for the three-dimensional Ising model at 
the ordinary transition, as revealed by scrutinizing more 
closely the early time behaviour. These deviations are of 
course due to the fact that for t — s <C s, i.e. outside of 
the dynamical scaling regime, the surface autoresponse 
(and therefore also the thermoremanent magnetization) 
exhibits a stationary behaviour. However, even in these 
cases the analytical curve again nicely match the numer- 
ical data for t — s > s. 

The situation at the surface is therefore similar to that 
encountered inside the bulktii in both cases no system- 
atic deviations from the theoretical prediction coming 
from the theory of local scale invariance are observed 
in the aging regime when analyzing the numerical data. 
This is a remarkable result, as it indicates that the scal- 
ing function of the autoresponse is completely fixed when 
knowing the values of two critical exponents: the scal- 
ing dimension x\ of the local order parameter [xi be- 
ing equal to x resp. x\ inside the bulk resp. close to 
the surface) and the scaling dimension Xi of the initial 
magnetization. At this stage one has to remember, how- 
ever, that renormalization group calculations^ yield in 
the bulk system a correction to the predicted behaviour 
(J3J at two loops. Similar corrections are then also ex- 
pected in a field-theoretical treatment of the surface re- 
sponse function in semi-infinite systems. It is an open 
problem why these corrections are not directly revealed 
when comparing the numerical data with the local scale 
invariance prediction. An obvious explanation would of 
course be that these corrections are in fact very small 
and therefore difficult to observe. There may, however, 
also be hidden a more fundamental problem. Indeed, ex- 
act results in one dimension yield the value Xc = 1 for 
the Glauber-Ising model at the critical point T = 0^ 
whereas one obtains from the time-dependent Ginzburg- 
Landau equation^ which is usually thought to describe 
the same system, the value Ac = 0.6006. This indicates 
that the one-dimensional Glauber-Ising model and the 
time-dependent Ginzburg-Landau equation belong to dif- 
ferent universality classes, which leaves open the possi- 
bility that similar problems arise in higher dimensions. 52 
Clearly, further investigations are called for in order to 
clarify this important issue. 

Before coming to the surface fluctuation-dissipation ra- 
tio in the next Section, let me briefly comment on the 
range of waiting times accessed in the present study. In- 
deed, a reader familiar with the investigations of aging 
phenomena in spin glasses may be surprized by the ap- 
parently small values of s retained here. In critical ferro- 
magnets the decay of, for example, the thermoremanent 
magnetization is very fast, due to the large values of Xr/z 
(being for the three-dimensional Ising model 1.73 in the 
bulk and 2.10 at the surface), whereas in spin glasses the 
decay is much slower with Xr/z typically of the order of 
0.1 — 0.4. This slow decrease makes it possible to access 
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much longer waiting times in spin glasses. On the other 
hand, however, the dynamical scaling regime is entered 
for spin glasses only for large values of s, whereas in ferro- 
magnets dynamical scaling behaviour is already observed 
for rather small values of s, as illustrated for instance in 
the different figures discussed in this work. 



V. SURFACE FLUCTUATION-DISSIPATION 
RATIOS 

One of the central quantity in the study of systems 
far from equilibrium is the fluctuation-dissipation ratio 
X(t, s) given by Eq. 10 £ In the context of critical sys- 
tems one is especially interested in the limiting value 



X°° = lim ( lim X(t, s] 

S >oo \t >oo 



(24) 



which characterizes the different universality classesiS^ 
For example, one obtains in the case of Glauber dy- 
namics the value X°° = 0.33 for the two-dimensional 
Ising model^^i whereas in three dimensions one has 
X°° w 0.4£ In the exactly solvable one-dimensional case 
this value is X°° = 1/2 at T = Oj&iS when starting from a 
uncorrelated initial stated This limiting value can also 
be obtained from the thermoremanent magnetization. 9 
Indeed, in the scaling limit the scaling functions fc 
and fu of the autocorrelation and response functions 
only depend on t/s, yielding M(t,s) = M(C(t,s)) and 
therefore^ 



X c 



lim 

c — >0 



T C M(C) 
hC ' 



(25) 



Fluctuation-dissipation ratios have up to now only been 
studied in bulk systems. 

We define the surface fluctuation-dissipation ratio in 
complete analogy to the bulk case, but with surface quan- 
tities replacing the bulk ones: 



X 1 (t,s)=T c R 1 (t,s)/ 



d s 



(26) 



The limiting value of the surface fluctuation-dissipation 
ratio is then given by the expression 



lim 

Ci — >o 



hCi ' 



(27) 



where C\ and M\ are the surface autocorrelation function 
and the surface thermoremanent magnetization studied 
in the preceding Sections. 

Figures |Ht and[SjD display the parametric plot of the 
rescaled thermoremanent magnetization s 2xi ^ z T c Mi/h 
versus the rescaled autocorrelation function s 2xi ^ z C\ 
for two different cases: the pure two-dimensional 
semi-infinite Ising model and the semi-infinite three- 
dimensional Ising model at the special transition point. 
In both cases the data points for different waiting times 
show the expected perfect scaling. A similar good 
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FIG. 8: The rescaled surface thermoremanent magnetiza- 
tion versus the rescaled autocorrelation function: (a) perfect 
two-dimensional semi-infinite Ising model with A — and (b) 
the three-dimensional semi-infinite Ising model at the special 
transition point. The dashed lines indicates the quasiequilib- 
rium regime at early times, whereas the slope of the full lines 
is given by the limiting value Xf° of the surface fluctuation- 
dissipation ratio, see Eq. 1251 and Table ITTT1 



scaling is observed for all studied systems in two and 
three dimensions. At early times, t — s <C s, a quasi- 
equilibrium behaviour, where the fluctuation-dissipation 
theorem holds, is indeed observed, as indicated by the 
dashed lines with slope 1. In the limit C\ — ► the 
slope of these parametric plots yield the asymptotic val- 
ues Xf° of the surface fluctuation-dissipation ratios, as 
indicated by the full lines. Alternatively, one may also 
investigate X\ as function of s/t, which yields the limit 
value X^° for s/t — > 0. Both approaches give com- 
parable results. The limiting values determined for the 
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TABLE III: The limiting value of the fluctuation- 

dissipation ratio obtained from Eq. 12511 for the different mod- 
els. HvL: Hilhorst-van Leeuwen model, OT: ordinary transi- 
tion, SP: special transition point. 



HvL 


A 




AT 




0.50 


1/4 


0.31(1) 




0.25 


3/8 


0.33(1) 







1/2 


0.37(1) 




-0.25 


5/8 


0.40(2) 




-0.50 


3/4 


0.43(2) 


d = 3 


OT 


1.26 


0.59(2) 


d = 3 


SP 


0.376 


0.44(2) 



different cases are gathered in Table ITTTl In the case of 
the Hilhorst-van Leeuwen model the quality of the exist- 
ing data do not allow a reliable determination of Af° for 
xi > 3/4. The main message of Table ITTTl is that Xf° de- 
pends on the surface scaling dimension: for the Hilhorst- 
van Leeuwen model a variation in the strength of the 
defect amplitude A changes the value of Af°, whereas 
in the three-dimensional Ising model the different uni- 
versality classes yield a different value for this limiting 
value. Assigning an effective local temperature via the 
equation T e ff = T c /X^°, one has to remark that (a) the 
surface effective temperature in general differs from the 
effective temperature T c /X°° inside the bulk, and (b) 
the surface effective temperature depends on the surface 
quantities (as for example the strength of the surface cou- 
plings) which determine the different surface universality 
classes. 



VI. CONCLUSION AND OUTLOOK 

The aim of the present work has been to extend the 
analysis of aging processes in critical systems far from 
equilibrium to the semi-infinite geometry, paying special 
attention to the behaviour of surface quantities in the dy- 
namical scaling regime. As shown in the preceding Sec- 
tions, surface autocorrelation and autoresponse functions 
display in this aging regime a scaling behaviour similar 
to that observed in the bulk, but with exponents differ- 
ing from those encountered inside the bulk. This scaling 
behaviour has been investigated numerically in the two- 
and three-dimensional semi-infinite Ising models as well 
as in the Hilhorst-van Leeuwen model. The latter model 
is a semi-infinite two-dimensional Ising model with an 
extended surface defect which has the property that the 
surface scaling dimension varies continuously as a func- 
tion of the defect amplitude. It is this property that has 
made possible a systematic study of the effects of surfaces 
on aging processes. 

It is remarkable that the scaling functions of the sur- 
face thermoremanent magnetization are found to agree 
with the analytical predictions coming from the theory 
of local scale invariance. 14 This extends the range of ap- 



plicability of this theory to semi-infinite systems. At this 
stage field-theoretical computations of scaling functions 
are called for. Remembering the situation in the bulk sys- 
tems, one may expect that renormalization group calcu- 
lations also yield corrections to the predicted behaviour 
(|22|1 in semi-infinite systems. This (expected) discrep- 
ancy between the observed agreement of the numerical 
data with the theoretical prediction on the one hand and 
the appearance of correction terms in a field-theoretical 
treatment on the other hand is an open problem which 
warrants more attention in the future. 

Similar to what has been done in the bulk systems^ 
field theoretical calculations should also permit to com- 
pute the limiting value Xf° of the surface fluctuation- 
dissipation ratio. The dependence of X%° on the surface 
scaling dimension, as found in this work, again illustrates 
that the limiting value of the fluctuation-dissipation ratio 
permits to characterize the different (surface) universal- 
ity classes. 

Possible extensions of the present work include the 
study of nonequilibrium semi-infinite systems in cases 
where the order parameter is conserved (Kawasaki 
dynamics)^ One would then be dealing with the semi- 
infinite extension of model B (in the classification of Ho- 
henberg, Halperin and Ma) , whereas in this work only the 
semi-infinite extension of model A has been considered. 
Indeed, the study of semi-infinite critical systems evolv- 
ing with Kawasaki dynamics gives rise to new questions. 
It has been shown that in the semi-infinite extension of 
model B with a conserved bulk order parameter the pres- 
ence of nonconservative surface terms, leading to a non- 
conserved local order parameter in the vicinity of the sur- 
face, yields a different dynamic critical universality class 
as compared to the case where these terms are absent. 
These two universality classes share the same critical ex- 
ponents but are characterised by different scaling func- 
tions of equilibrium dynamic surface susceptibilities^^ 
It is therefore tempting to ask whether this different 
surface universality classes also show up when studying 
nonequilibrium quantities. Especially, it would be in- 
teresting to investigate the scaling functions of the au- 
toresponse in these cases and to compare them with the 
prediction (|22|l . A second obvious extension would be 
the study of semi-infinite systems at temperatures below 
T c where phase-ordering takes place. Dynamical scaling 
behaviour is also observed in bulk systems in the low tem- 
perature phase and a similar behaviour is also expected 
for semi-infinite systems. Work along the sketched lines 
is planed for the future. 
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(t 2 d t +s 2 d s + 2&i + 2C 2 s) Ri (t, s) = (A.2) 



APPENDIX 



with the solution 



In this Appendix I derive the expression for the 
scaling function of the surface autoresponse under the 
assumption that local scale invariance 14 holds. I hereby 
closely follow the derivation of the scaling function of the 
bulk autoresponse given in Ref. Q. 

Consider the local two-point function TZ±(t, s; p\, pj) 
which describes the response of the system at time t at 
the surface site p\ to a perturbation acting at time s at 
the surface site f>2- As we are considering aging systems, 
time translation invariance is broken. However we still 
have space translation invariance along the surface, and 
therefore 1Zi(t, s; p\, pi) = 1Zx(t, s; p) with p = p\— pi- 
We furthermore require that TZ,\ transforms covariantly 
under scale (Xq) and special conformal (Xi) transfor- 
mations, i.e. Xo IZi = X\ IZi = 0, where the differen- 
tial operators Xq and X\ are explicitly given in Ref. Q. 
The differential equations for the surface autoresponse 
Ri(t, s) = Hi(t, s; 0) are obtained by setting p = \p\ = 0, 
yielding 



C2-C1 

R 1 (t,s)=r [- s ) (t 0(t- a). (A3) 



Here, £1 and £2 are two exponents left undetermined by 
the theory, whereas ro is a normalization constant and 
O the step function. The exponents Ci and C2 are fixed 
by comparing expression (|A.3Jl with the expected scaling 
behaviour 1)13(1 and (|14[1 . leading to the final result: 



• ,v l+2xi/»-Ai/« 

./?,(/..) = ,•„ ( - \ ( t _ s )-i-2^ e( i_ s) . 

(A.4) 

This can be written in the form 



R 1 (t,s) = s- 1 - 2x ^f Rl (t/s) 



(A.5) 



{td t +sd s +(l+( 2 )Rl(t,s) 







(Al) 



where the scaling function f Rl (y) is given by the expres- 
sion 
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